Preprint numbers: ITP-UU-12/32, SPIN-12/30 



On antiscreening in perturbative quantum gravity and resolving the Newtonian 

singularity 

Anja Marunovic a * and Tomislav ProkopecQ 
a University of Zagreb, Faculty of Electrical Engineering and Computing, 
Physics Department, Unska 3, 10000 Zagreb, Croatia 
b Institute for Theoretical Physics and Spinoza Institute, 
Utrecht University, Leuvenlaan 4, 3584 CE Utrecht, The Netherlands 

We calculate the quantum corrections to the Newtonian potential induced by a massless, non- 
minimally coupled scalar field on Minkowski background. We make use of the graviton vacuum 
polarization calculated in our previous work and solve the equation of motion non-perturbatively. 
When written as the quantum-corrected gauge invariant Bardeen potentials, our results show that 
quantum effects generically antiscreen the Newtonian singularity oc 1/r. This result supports the 
| point of view that gravity on (super- )Planckian scales is an asymptotically safe theory. In addition, 

we show that, in the presence of quantum fluctuations of a massless, (non)minimally coupled scalar 
P-t field, dynamical gravitons propagate superluminally. The effect is, however, unobservably small and 

1 it is hence of academic interest only. 
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I. INTRODUCTION 



In Ref. [1| we considered the graviton vacuum polarization induced by the one-loop quantum fluctuations of a 
qj _ nonminimally coupled, massless scalar field on Minkowski background in the Schwinger-Keldysh formalism. This 
i <~| allowed us to compute the large distance, leading order, quantum correction to the scalar gravitational (Bardeen) 
l— — potentials, which in longitudinal gauge reduce to the Newtonian potentials. We used the Schwinger-Keldysh formalism 
I . in direct (coordinate) space which allowed us to study time transients associated with switching on scalar vacuum 
1 fluctuations at a definite moment in time. Such transients can be induced, for example, by symmetry restoration 
Q\ \ that renders the scalar field massless. Moreover, just as in the presence of classical matter, the two quantum Bardeen 
f — potentials differ when the effects of quantum (scalar) matter fluctuations are included. 

In Ref. [l| we adopted the effective approach to quantum gravity (see e.g. Refs. 0, H|), according to which quantum 
gravitational effects can be reliably calculated by perturbative methods, provided one focuses on low energy and 
Q\ * momentum phenomena. 1 By now, there is quite a long history and extensive literature on long-range, quantum 
corrections to the Newtonian potentials. Such studies were pioneered by Donoghue @, and followed by many 
others [1, with differing results. The problem was again revived in the 2000s in Refs. |6l-fl2l|. where also some 
of the graviton vertex corrections were included. The Schwinger-Keldysh formalism was for the first time used in 
i> Ref. [13J . Related older works include [14l . |15| where the one-loop graviton vacuum polarization induced by photons 
J^j [ on a Minkowski background was calculated. 

Here we extend our work in [1] and present a resummation scheme by which we consider short and large scale 
screening effects of gravitational potentials induced by quantum fluctuations of a massless nonminimally coupled 
scalar. More precisely, we make use of the one-particle irreducible (1PI) resummation technique and arrive at the 
result that the Newtonian potentials exhibit a strong antiscreening at length scales of the order of the Planck scale. 
Surprisingly, there are also long range oscillatory quantum corrections to the Newtonian potentials. Even though 
these corrections are present at large distances, they are high energy corrections, since oscillations occur at the Planck 
length scale. 

Admittedly, these results do not belong to the realm of effective quantum gravity, questioning their reliability. 
Still, we feel compelled to present them, and believe that they deserve some attention. Indeed, this work is not the 
first attempt to use a resummation technique to investigate short scale effects in quantum gravity. The main other 
contender is known as the effective action approach to quantum gravity [16[ , and it is usually implemented within the 
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1 By low energy/momentum phenomena we mean the phenomena whose energies E and momenta p = \\p\\ are well below the Planck 
scale, E <C Ep, p <C Ep/c, where Ep = mp c 2 = yfhc~>/G N ~ 1.2 X 10 19 GeV, where mp is the Planck mass, c is the speed of light, 
Gjv denotes the Newton constant and h = 1.055 X 10~ 34 Js is the reduced Planck constant. 
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asymptotic safety program for quantum gravity, and we shall now recall the basic idea behind this approach. 

The asymptotic safety program in quantum gravity was initiated by S. Weinberg [17| . where he argued that a 
weaker requirement than the usual renormalizability may be imposed on quantum gravity. Recently, this program 
has attracted a considerable attention among researchers in quantum gravity pH . Il9j . If quantum gravity were an 
asymptotically safe theory, it would exhibit a few dimensional critical hypersurface characterizing the running of 
relevant couplings (that are associated with relevant operators). The running of these couplings 'emanate' from an 
ultraviolet (UV) fixed point, at which their value has a definite value that can be, in principle, determined. The claim 
is that the set of relevant operators consists of unity (the corresponding coupling being the cosmological term), the 
Ricci scalar (the corresponding coupling being the Newton constant) and possibly some quadratic curvature invariants. 
There is a limited evidence that all other (higher derivative) operators are irrelevant, in the sense that the UV fixed 
point 'repels' the corresponding couplings as one runs towards the ultraviolet. The value of these irrelevant couplings 
is (at least in principle) uniquely determined by the requirement that the physical theory runs towards the UV fixed 
point, as the running scale approaches infinity. Thus, an asymptotically safe quantum gravity consists of several 
(relevant) couplings whose value ought to be determined by measurements, for example in the far infrared, while all 
other coupling parameters are given by the requirement of asymptotic safety. Within this framework it is not true 
anymore that there are infinitely many counter-terms appearing in perturbation theory, whose coefficients (couplings) 
must be determined by (an infinite number of) measurements. In this sense the asymptotic safety program resolves 
the problem of non-renormalizability of quantum gravity. What remains to be shown is that quantum gravity is 
indeed an asymptotically safe theory, which is far from a trivial pursuit. 

Indeed, the proof that quantum gravity is asymptotically safe would require solving quantum gravity non- 
perturbatively, which is very hard. Instead, researchers have focused on developing perturbative and non-perturbative 
methods (for reviews, see [IH, UH), hoping to collect enough evidence in support of the asymptotic safety hypothesis. 
The original evidence due to S. Weinberg [17] is based on an e-expansion around 2 spacetime dimensions. Here, e = 2 
is not a small number in four spacetime dimensions, making the whole approach questionable. Second approach is 
based on functional renormalization group equation for the gravitational effective action pioneered by Reuter 
While the starting equation is an exact (functional) differential equation for the effective gravitational action, the 
actual approximation scheme is limited to an a priory chosen set of local operators, truncated at some order (usually 
fourth) in space and time derivatives. The criticism of the method is twofold. Firstly, the operators are chosen a 
priory based on a derivative expansion which is essentially a low energy expansion, while the application is extended 
to the ultraviolet sector. Secondly, the operators of canonical dimension four or higher generically violate unitarity, 
questioning thus the validity of the whole (non-perturbative) scheme. Yet, the results (if true) are fascinating: they 
suggest existence of an ultraviolet (UV) non-Gaussian, non-trivial, fixed point in quantum gravity. Close to this UV 
fixed point the Newton constant runs with a scale /i as Gn{^) — 9*/fJ- 2 (as /i — > oo) , while the cosmological constant 
runs as A(/i) ~ A*^i 2 , where g* ~ 0.27 and A* ~ 0.35 are calculable (non-vanishing) real constants. 

In this work we present a different resummation scheme, such that antiscreening occurs as a result of resumming the 
one-loop 1PI diagrams for the gravitational potentials. In contrast to the local operator approach mentioned above, 
where a derivative expansion is used, our resummation is based on a non-local correction to the graviton vacuum 
polarization, and thus does not suffer from the above problems. The hope is that, if we choose the coefficients of the 
finite counterterms to be zero, our method will not violate unitarity. Furthermore, our method does not violate any 
of the symmetries of the theory. Namely, we linearize in the gravitational field, and in our analysis we maintain the 
linear version of diffeomorphism invariance (also known as gauge invariance). 

The paper is organized as follows: in section II we present a self-consistent solution to the one-loop 1PI effective field 
equations for the gravitational potentials. In section III we show how the lowest energy (circular) orbits of classical 
point particles get affected by the long range potential that we found. In section IV we discuss our results and give 
an outlook for future work. 

Unless stated explicitly, we work in natural units where h = 1 = c. 

II. RESUMMATION 

The quantum corrected equation of motion for the linearized Einstein equation (see Eq. (39) in l]) is of the form: 

W^W + / dV Un£] (x; x')h^(x') + 0((V) 2 ) = y^MJ 3 ^), (1) 

where [/u,IE^] (x; x') is the retarded graviton vacuum polarization tensor of the graviton field perturbation h^(t, f), 
M is a (static) particle's mass and k 2 = IQttGn, Gn is Newton's constant. A massless, nonminimally coupled scalar 



+ c^^^ + 




FIG. 1: The resummed one-loop contribution from scalar vacuum fluctuations to the graviton propagator. 
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FIG. 2: The resummed contribution from the non-local (cubic) vertices to the graviton propagator. The resummed loop can 
be written as a sum of the one- loop, two-loop, etc., contributions. 

field on Minkowski background induces at one-loop the following renormalized graviton vacuum polarization tensor, 
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where Ar — \\x — x'\\, At = t — t' (x, x', t and t' denote spatial and time coordinates, respectively), £ measures the 
coupling strength of the scalar filed ip to gravity through the Ricci scalar R, 2 and Lp Vpa denotes the Lichnerowicz 
operator on Minkowski background: 3 



with the explicit form of the operators: 
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d{ P Va)(pd v) - ^iip(pV*)ud 2 ~ ^>l P <ydpd u , 

q2 _ r jjSQ^Q s an( j indices in (p} are lifted with the Minkowski metric, rf" = diag(— 1, 1, 1, 1). 

In this paper we solve the 1PI equation of motion for the graviton one-point function (JlJ exactly. This corresponds 
to resumming all 1PI one-loop diagrams, illustrated in figure [TJ where we show all: the non-local contribution arising 
from cubic graviton-scalar vertices, the local contribution from quartic scalar-graviton vertices, and the counterterm 
contribution. In order to illustrate how the resummation works, in figure[2]we show the resummed non-local diagrams 
(arising from cubic vertices). In principle, the contribution from the local diagrams needs to be resummed as well. But, 
in a flat spacetime background they contribute only as power-law divergences, their contribution gets automatically 
subtracted in dimensional regularization, and thus they do not contribute 

To proceed, note first that it helps to insert instead of the source on the right hand side of (JXJ) the classical equation 
of motion 

„2 
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The covariant action of the nonminimally coupled, massless scalar field ifi is: 



: Ei 



2 , rr,, -r,, r 

3 Recall that, up to boundary terms, the quadratic Einstein-Hilbert action can be written in terms of the Lichnerowicz operator as 
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where superscript c stands for classical. This enables us to completely remove the Lichnerowicz operator from Eq. ([T]) 
and thus relate the quantum gravitational field to the classical potential through an integral equation: 



h a ^ c \x) = h a0 {x) 



n a 

7680^ 1 p ° 



pa 



J d 4 x'Q{At 2 - Ar 2 )9(Ai) [ln(^ 2 (Ai 2 - Ar 2 )) - l] h pa {x') 



(6) 



This form of the equation suggests that the vacuum polarization can be thought of as a (non-local, spacetime 
dependent, causal) wave function renormalization of the graviton h a ^(x). In the Appendix we show how to 
solve this equation. The general solutions in momentum space are given in Eqs. (|80H81[) . Here we are inter- 
ested in a static response to a point mass, for which the classical Bardeen potentials in momentum space are, 
$(c) = ,j/( c ) = -^TrGjvAf/fc 2 ] x 27T(5(fc°), k = \\k\\ (see Ref. [l]). Hence, a static response to a point mass M is 
obtained simply by taking the limit k° — > of Eqs. 
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Since for the vector and tensor perturbations classical fields vanish, Eqs. ([7411751) imply that, at one-loop order, also 
the quantum corrected fields vanish, hf = = hf^ ' . 

Upon performing an inverse Fourier transformation and setting x = klpi, jl = ^,lpie~ lE , y = r/lpi (where Ipi = 
■\/G~n is the Planck length), we arrive at the solution for the gauge invariant scalars in position space: 
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As such these integrals cannot be computed analytically. However, it is feasible to examine their behaviour for small 
r and for large r. 

For small r and for nonconformal coupling (£ ^ 1/6) both integrals are proportional to r with the constant of 
proportionality depending on the coupling strength £ only. This is the main result of this section as it clearly shows 
that the 1/r singularity is screened on Planck scales by quantum effects (see figures [3] and S] for rs/lpi = 10 and 
figures [5] and O for rs/lpi = 0.05, where r$ — 2GjyM is the Schwarzschild radius). Nevertheless, in the case of the 
conformal coupling (£ = 1/6) the second integral in the above expressions reduces to a constant while the first one is 
proportional to r. This means that for the conformal coupling there is no screening of the Newtonian singularity (in 
this case quantum effects do not increase the order of singularity at the origin), which means that quantum effects 
are the weakest in the conformally coupled case. 

For large r, the first integral, denoted by I\{y) 



h(y) 



dx sin(xy) 



1 



1 + x In 

1 ' 1207T 
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4 In general, the classical solution h a P( c ' in Eq. J6j includes also gravitational waves, which are solutions to the homogeneous equation, 
L^i/pahf"^ 3 ^ (x) = 0. They are not of main interest here, since we are mainly interested in the gravitational response to a static 
point-like massive particle. 
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is rather simple and yields a constant I\{y) = tt/2, while the second one, we denote by hiy), 

f °° sm(xy)dx 1 

h(y) = / ; 7—. v ( 13 ) 

Jo X 1 _ (1 _ 6e) 2^. ln ^ e 2 7B 

is more delicate and needs to be treated heedfully. As £ ^ 1/6, the denominator in this integral has poles on the 
real axis denoted as ±x c . In order to calculate principal part of the integral, we approximate the log-function by a 
constant evaluated at x c . This is justified by the fact that the log-function is the largest (but finite) at the pole. In 
this approximation the integral becomes 

lf°° ( 1 1 2\ 
h{y) = —z / dxsm(xy) + — , (14) 

Z Jo \X X c X + X c X J 

where x c — x c (£,fj,) is the (positive, real) solution of the transcendental equation 

The approximate integral (| 14[) can be expressed in terms of the sine integral Si(x) and the cosine integral Ci(x), 
defined by: 

= _ f°° sMz)dz = f°° cos( Z )dz 



Using asymptotic expansions of Si(x) and Ci(a;), we get the following approximate expression for J 2 in the limit of 
large r: 

h{y) ~ | [1 - cos(at c y)] . 

This result needs to be checked against the numerical one. It turns out that one gets agreement, provided one modifies 
the coefficient in front of the cosine function by a function that weakly depends on £. So a better approximation for 
asymptotically large r is: 

7 2 (y) = |[l-a(0cos(x c y)] , (17) 

where a(£) ~ 1. When this and I\ = tt/2 are inserted into ([T0l - [TT|) one arrives at the following asymptotic form for 
the Bardeen potentials: 



{l + ia(e)cos( Kc r)} (18) 
^^{l-ia(£)cos( Kc r)}, (19) 



where yx c — n c r, k c — x c Jlpi. Notice first that, unlike in the classical case, the one-loop scalar quantum fluctuations 
generate different Bardeen potentials at large r. Second, Bardeen potentials exhibit oscillatory behaviour at large r. 
These oscillations occur at the Planck scale, and their amplitude is about 1 /3 of the classical potential (modulated 
by a weakly dependent function of £ of order 1). Eqs. (TT8"|) and (TlT)!) comprise the second main result of this section. 

In figures [3] - [6] we show the numerically evaluated quantum corrected Bardeed potentials as a function of the 
distance from the point mass M. While figure|3]and figure@]show the Bardeen potentials for a heavy mass, M ^> mpi, 
figure [5] and figure |6] show the Bardeen potentials induced by a light mass, M -C rripi (which holds for the mass of 
all known elementary particles). The figures show that, when quantum fluctuations are included (solid lines), the 
Newtonian potentials (short-dashed lines) get antiscreened such that the Newtonian singularity oc 1/r gets resolved. 
Furthermore, for large r, there is very good agreement between our semi- analytical expressions for the Bardeen 
potentials flTBl and (|19p and the curves obtained by numerical integration. Next, we show how the potentials depend 
on the rs/lpi parameter. In general, if the Schwarzschild radius is greater than the Planck length then an event 
horizon forms (a property characterizing black holes), in which case the Bardeen potential *f> crosses —1/2. This is 
exactly what we see in figure 2J However, if the Schwarzschild radius is less than the Planck length (as it is the case 
with all known elementary particles) then no event horizon forms, i.e. ^> > — 1/2 for all r. This feature is clearly 
manifested in figures [S] and [SJ 
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FIG. 3: The quantum corrected Bardeen potential <&, obtained numerically, is represented by the solid line, long-dashed line 
represents its analytical expression (|18[) for large r behaviour, short-dashed line shows Newtonian potential. Values of the 
parameters are: £ = (or £ = 1/3), p. = 1 and rs/lpi = 10. 




FIG. 4: The quantum corrected Bardeen potential ^, obtained numerically, is represented by the solid line, long-dashed line 
represents its analytical expression (|19[) for large r behaviour, short-dashed line shows Newtonian potential. Values of the 
parameters are: £ = (or £ = 1/3), p. = 1 and rs/lpi = 10. 



We have seen that quantum fluctuations antiscreen gravitational potentials in the sense that they reach a finite value 
when r — > 0. For small r the potentials are linear in r, i.e. $ ~ <£>o + &Q r , & — + ^o r - The resulting gravitational 
force, F = — V$ = — <&' r/r is constant and negative, implying that the force reaches a maximum amount. This is 
to be contrasted with the Newton force, which grows without a limit as oc 1/r 2 at small r. Let us now consider the 
Riemann curvature tensor R paliV . When linearized in h^v, 
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Taking account of symmetries of the Riemann tensor, the nonvanishing components of (|20[) can be written in terms 
of gauge invariant components of h^ v only as, 
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FIG. 5: The quantum corrected Bardeen potential 3>, obtained numerically, is represented by the solid line, long-dashed line 
represents its analytical expression (|18p for large r behaviour, short-dashed line shows Newtonian potential. Values of the 
parameters are: £ = (or £ = 1/3), jx = 1 and rs/lpi = 0.05. 
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FIG. 6: The quantum corrected Bardeen potential obtained numerically, is represented by the solid line, long-dashed line 
represents its analytical expression (|19p for large r behaviour, short-dashed line shows Newtonian potential. Values of the 
parameters are: £ = (or £ = 1/3), jl = 1 and rs/lpi = 0.05. 

When applied to the problem at hand (hf? = 0, = 0, nf = 0), these reduce to, 

— ~)^'o ' Roiji = o , 

R ljkl = 4(^-2lfe)*' 0! (22) 
from which we can obtain the Ricci curvature tensor and scalar, 

R 00 = -%, ~\ L ) '''•> ^Wo. fl = --$' -V . (23) 

r V r r d / V r r 6 / r r 

From these equations we see that a curvature singularity remains. It is however reduced from the standard curvature 
singularity of a black hole, R p(Tl j,vR pa>lv = (48G^,A/ 2 )/r 6 (which is present only in the Weyl part of the Riemann 
tensor since = 0) to much milder curvature-squared singularities, 

= 4( 3$ o 2 + W &o + n*' Q 2 ), R pa ^R^ = ^(<P' 2 + 3% 2 ), (24) 
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which for small r are of the order ~ 1/r 2 . It would be of interest to check whether this behavior persists at two 
loop order, or perhaps the screening of gravitational potentials is more efficient when higher loop effects are taken 
into account. Such an analysis would increase our understanding of the strength of gravity at very small (Planckian) 
scales. 



III. PLANCK SCALE MODULATION OF CLASSICAL ORBITS 



In this section we analyze how the quantum fluctuations stored in the Bardeen potentials (TT5)) and (|19p affect 
circular motion of a classical particle. Hence we start with the geodesic equation 

d 2 x» u dx a dx 13 , x 

For this section we shall work in the longitudinal (Newton) gauge, in which the Bardeen potentials reduce to the 

(longitudinal) potentials, $ —> <fi and \& — > ip. The metric written in terms of the these potentials (see Eqs. (fT8|) 
and ([II])) is: 

ds 2 = -(1 + 2(f))dt 2 + (1 - 2i}))dr 2 + r 2 (d6 2 + sin 2 8dtp 2 ). (26) 
From the following two components of the geodesic equation 

d 2 t d r 4> dr dt 

d}? + T+24>dXdX = °' ( 7) 
d 2 ip 1 „ v,„ , s 9i dtp dr „ cos 9 dO dip „ 

W + W^W) dT [(1 - 2 ^ )r ] ^Tx + 2 ^oTxl = °' (28) 

we extract the two Killing vectors: 

= (-(1 + 20), 0,0,0), (29) 

Ffj, = (O,O,O,r 2 (l-2^)sin 2 0), (30) 

which lead to the conserved quantities (see e.g. (20|): 

'-'•l-f^ 1 E = ' (31) 

L = F M — =r 2 (l -2i/;)sin 2 0-^, i.e. -77^ = 0. (32) 



Here 15 is the conserved energy (per unit mass) and L is the conserved angular momentum. Since the direction of 
the angular momentum is conserved, without loss of generality we can set 6 = tt/2. From the velocity normalization 
condition u^u^ = — 1 and making use of the conserved quantities E and L we obtain an expression for the radial 
velocity squared: 



drV if E 2 L 2 



dXJ l-2r/>ll + 20 r 2 (l-2^) 
If we rewrite this equation in a slightly different form 



-1 . (33) 



1 /dr x 1 



2V,/A' - ]:iiir) 
where e — E 2 /2, the effective potential V^s (r) can be read off: 

L 2 1 _E 2 / 1 \ 

VcS W = 2r 2 (l-2^) 2 + 2(1-2^) + ~ V (1 + 20)(1 - 2^) J ' (35) 

Since we are interested in the large r behavior, we linearize this potential around small 4>, tp yielding: 

V eS (r) « E\cf -iP) + Ul + ^)+ip(l + ^f). (36) 
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The radii of the circular orbits are located in the minimum of this potential, which solves the equation: 

—T + o °(0 1 cos(Kr) 2£ 2 - 1 T )+ (r«) sin(«x) 2£ 2 - 1 =- } = 0. (37) 




GjqMr r 2 3 I \ r 2 / \ r 2 

The radii of the classical orbits are located at r c = GjyM/L 2 . This is what we get from the first three terms in 
the above equation in the large r limit (as the term oc 1/r 2 can be neglected). Second part of the above equation, 
multiplied by a(£), for r s = 2GnM <C L Cr, can be compactly written as 



as S is much smaller than nr (fix ~ r/lp\): 



sm{nr + 5) = 0, (38) 
2E 2 ~ 1 - 6L 2 /r 2 1 Zpi 



(39) 

r 



nr(2E 2 - 1 - 2L 2 /r 2 ) Kr 
From the 'quantization' condition (|38[) we obtain: 

K,r n +6 = nir n=l,2, ■■■ (40) 



If we now neglect S, from (|40|) it ultimately follows that the correction to the radii of the classical circular orbits due 
to the quantum fluctuations is: 

r n — — nir n = 1, 2, • • • (41) 

An interesting question is whether this phenomenon of a Planck scale 'quantization' can be observed. Here we made 
a very crude point particle approximation, which effectively treats (low energy) particles as classical point particles. 
This means that in our treatment point particles can be localized to a precision comparable or better than the Planck 
scale, which is obviously at odds with quantum theory. Indeed, the Heisenberg uncertainty relation tells us that a 
particle of a momentum p cannot be better located than Ax ~ h/p. This then means that, for particles with momenta 
and energies p,E <g; Ep\, the Planck scale orbits considered in this sections will appear smeared. To find out precisely 
what effect the Planck scale oscillations have on the motion of a quantum particle would require a complete quantum 
treatment of planetary orbits, which is beyond the scope of this paper. Of course, it is a very interesting question 
whether these Planck scale (but large distance) oscillations could be observed, and what precisely would be their 
experimental signatures. 



IV. QUANTUM GRAVITATIONAL RADIATION 

Here we discuss how the one-loop vacuum polarization from (nonminimally coupled massless) scalars ^ affects 
the propagation of gravitational waves. For simplicity, we perform here a perturbative analysis of the 1PI resummed 
equation ()75[) . and we shall work in direct space. This section is inspired by the recent paper of Leonard and 
Woodard [2l|, where it was shown that the photon one- loop vacuum polarization induced by the graviton vacuum 
fluctuations affects propagation of the photons by deforming its light cone. Here we shall see that the analogous effect 
is present in the case of dynamical gravitons immersed in a sea of vacuum fluctuations of matter fields. 

We begin our analysis by a simple quadrupole model, for which non- vanishing components of the conserved energy 
momentum tensor are: 

T Q0 (x,t) = 4-r(t)0i0,<y 8 (2) 

T 0i (x,t) = j tf e(t)a,-a 3 (f) 

T tJ {x,t) = l l3 5{t)8 3 {x), (42) 

where 1^ = J d 4 xe lk ' x Tij(x) is the transverse, traceless, momentum space moment of inertia, satisfying la = and 
dilij = and r(t) is a 'ramp' function dr{t)/dt = 0(f), and dQ(t)/dt — S(t). Of course, the stress energy tensor (j4*2"j) is 
an approximation to a realistic explosive event with a non- vanishing quadrupole, which has a finite spatial extension. 
The classical, gauge invariant equation for the (transverse-traceless) graviton is 

d 2 h\f TT (x) = -l&rdtTij (43) 
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and its solution for the stress energy tensor (142]) is the Einstein quadrupole formula, 

h\f TT (x) = 4Gjv4-fc^ = 8G JV 4-e(i)<5(x 2 ) , x 2 = -t 2 + r 2 , r= \\x\\ . (44) 



An inspection of Eq. (|T5[) shows that, at this order, quantum effects do not mix tensor perturbations with scalar or 
vector perturbations. Thus, Eqs. ([55H54]) can be rewritten for hf? as 



h$ TT (x) = hjftx) - T^\D iJU - 8 l3 D kl 



1 / 1\2 

6 +30 (^6) 



F M (x), (45) 



where 



F ki {x) = 2 J d 4 x'e{At 2 - Ar 2 )9(Ai) [ln(^ 2 (Ai 2 - Ar 2 )) - l] d' 2 h T J(x') , (46) 

where we have moved by partial integration two out of four derivatives to act on h^f . Eqs. (143114^)) can be solved 

perturbatively, by writing h\J — h kl ' + h kl , where h kl ' is the classical solution (l44l) and h kl ' is the quantum 
correction. Inserting this into (|45rfloT) results in, 



1 / 1\ 2 



Fu(x). (47) 



Inserting h k c ^ TT into (|46|) yields the leading order (one-loop) quantum correction to the graviton field. Making use 
of (g21), Eq. (gSJ) simplifies to, 

F kl (x) = -WTrGN^d 2 {Q(t 2 - r 2 )Q(t) [\n([i 2 (t 2 - r 2 )) - l] } . (48) 

Acting naively with the last two derivatives is problematic as it leads to singular expressions. This can be resolved 
by recalling that in the Schwinger-Keldysh formalism the retarded vacuum polarization originates from, 

[lnV4 + ) - 21n( M 2 4 + )] - [ln 2 (/i 2 4_) - 2\^ 2 x 2 + _)} \ 



Fki(x) = -16iTG N I l3 d 

32G N ~ (\n{[i 2 x 2 ++ ) ln^ 2 ^.) 
= hi i n o 



V 



Am 

2 x 2 ++ ) \d.(h 2 x 2 

2 l2~ 



(49) 



where x\ + — —(\t\ — it) 2 + r 2 , x 2 = — (t + ie) 2 + r 2 . As a final step we insert this result into Eq. (|47|) . Two 

important observations are in order. Firstly, the operator Dij defined in Eq. Q yields zero when it acts on a traceless, 
transverse tensor, implying that the term in (|45p containing the Ricci-scalar coupling £ does not contribute to the 
graviton propagation, i. e. a minimally and nonminimally coupled scalar contribute equally to graviton propagation. 
Secondly, the operator Dijkl defined in Eq. (U) yields a very simple result when it acts on a transverse, traceless 
tensor, DijkihJ^f = — (l/2)<9 2 /i^ T . With these observations, it is now easy to evaluate the leading quantum correction 
to the graviton field (HTl) . The result is, 

Mtt ( s_ 2G 2 ~ r 1 1 \_ 2G 2 N ~ d / 1 1 \_ 4G 2 N ' 



ft« (*) = ~T^^ W i-\ = TFT J «1TI \- —\ = -TF^^f2 (*) e (* - r ) • ( 5 °) 

J 15tt 2 i l_ x\ + x\_ J 157r 2 i dr 2 [ x z ++ x\_ J 15n dr 2 

Combining this with the classical graviton field (|44|) gives, 

h TT {x) = 8GnIij ^ _ e{t)5{t 2 _ r 2 } + = g^.^ - r 2 + ^) + 0(G%) . (51) 

Just as in Ref. [UJ, where the authors found that the photons in presence of graviton quantum fluctuations move 
slightly superluminally, we have found here that the gravitons in presence of scalar quantum fluctuations move also 
slightly superluminally 5 , i.e. they move on a shifted light cone given by, 

" 3 = t2 + S- < 52 » 



5 If one defines the graviton speed as the ratio r/t, then Eq. I|52|l implies that gravitons move slightly superluminally. However, if one 
looks at the local speed dr/dt = t/r (which is more in the spirit of group velocity), then the local speed is subluminal, and as t — > 
it even reaches zero. This is of course just a curiosity, since the lightcone deformation i r >'2\i becomes perturbatively trustable only for 
super-Planckian times, t 2> tpi. 
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FIG. 7: The classical and quantum corrected light cones of dynamical gravitons. 

Curiously, the effect docs not depend on £, i.e. it does not depend on whether the scalar is minimally or nonminimally 
coupled. This means that the effect (1521 comes entirely from the kinetic term cubic vertex —h^ v {d ll ip)(d v ip)/2. The 
classical light cones, as well as the quantum-shifted light cones, are shown in figure [Jj Just as in the case of the 
photons, this graviton effect is very tiny (in fact, it is 40 times smaller than the effect found for the photons in 21]), 
and thus it is uncertain whether it will ever be observed. From the structure of the light-cone in figure it is 
clear that the superluminality is not cumulative (it does not increase with distance). On the contrary, it becomes 
weaker as the distance increases. This can be understood as follows: the effect on larger scales is built primarily 
from quantum fluctuations on that scale, which have smaller amplitude than fluctuations on a smaller scale, and v. v. 
Yet, the effect is of some academic interest, as it is for the first time that a superluminal propagation of gravitons 
is inferred from quantum matter fluctuations. As a curiosity, we mention that a similar deformation of light cones 
has been observed in Hermitean Gravity of Ref. [22| , where a slight superluminal behavior was observed for particles 
with negative 4-momentum squared, — p 2 = \(E/c) 2 — p 2 ] > 0. Our effect is also different from the superluminal 
photon motion observed in curved geometries |23l425| . where the effects were calculated in a gradient expansion. No 
such approximation has been used here. Moreover, our results are manifestly gauge invariant, as they are calculated 
within a gauge invariant formalism. 

V. SUMMARY AND DISCUSSION 

In this paper we perform one-particle irreducible (1PI) resummation of the one-loop vacuum fluctuations of non- 
minimally coupled, massless, scalar matter, and study its (backreaction) effects on classical point particles moving on 
bound orbits. Our main result $TE\ indicates that gravity gets strongly antiscreened on (super-)Planckian energies, 
supporting the hypothesis of asymptotic freedom, according to which gravity becomes weak on super-Planckian scales. 
Furthermore, we study how vacuum fluctuations affect graviton propagation. We find that gravitons moving in a sea 
of vacuum fluctuations are sped up, resulting in a slightly superluminal motion (see figure UJ. However, the effect is 
too weak to be observable. 

Antiscreening on short scales is such that the Bardeen gravitational potentials remain finite everywhere, also 
arbitrarily close to the origin. However, the Riemann curvature tensor still exhibits singularity at the origin, albeit 
this singularity is much milder {B, iLvpa By ,vpa oc 1/r 2 ) than the usual singularity of black hole space times (for which 
RfivpaR^ vpcr °c l/r e ). In addition, the Ricci scalar and Ricci tensor are also singular, which is not surprising since 
vacuum fluctuations act as matter. Since for (elementary) particles of sub-Planckian mass the Bardeen potentials 
remain small everywhere (see figures [5] and [5]) , our results may be trustable when applied to known elementary 
particles. For heavy particles (whose mass is comparable to or larger than the Planck mass), we also find qualitatively 
the same strong antiscreening. However, our results cannot be trusted here since, at short distances (closer to the 
would-be event horizon), the amplitude of the gravitational potentials becomes of the order one or larger (see figures [3] 
and S]), invalidating our linearized approach to gravitational perturbations. So, in order to obtain more reliable 
results also in this case, it is required to extend our methods beyond linear order in the gravitational fields, or work 
in coordinates in which perturbations remain small everywhere. But, this task is left for a future work. 

A couple of remarks follows. It is believed that, just like Quantum Chromodynamics (QCD), quantum gravity is an 
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asymptotically free theory, in the sense that in the far ultraviolet it decouples from matter fields and becomes weak. 
There is some evidence for the asymptotic safety program that in the UV quantum gravity develops correlations as 
if it were effectively a two dimensional theory. This then implies that the corresponding spectrum is scale invariant, 
which may have relevance for generation of cosmological perturbations. Asymptotic freedom of QCD has lead to the 
conjecture (which has in the meantime been confirmed by numerical simulations) that, at high temperatures, the 
QCD matter is in a state of a quark-gluon plasma, in which most of the quarks (those of a sufficiently high energy) 
move essentially free. A quite exciting possibility is that a similar scenario may be realized with gravity. One can 
then imagine a collection of gravitationally weakly interacting very closely packed particles which would make dense 
objects which - due to the gravitational asymptotic freedom - would not collapse into black holes, or at least would 
screen curvature singularities at the origin of black holes. Whether this picture is consistent with our understanding 
of quantum gravity remains to be seen and requires further investigation. 

Our analysis in the Appendix contains results that are worth commenting on. The graviton vacuum polarization 
tensor of a massless nonminimally coupled scalar field has such a tensor structure that, when it acts on a graviton 
field, it yields gauge invariant field components. This result motivates the following conjecture: 

Gauge invariance dictates the structure of the graviton vacuum polarization tensor induced by any massless 
matter or gravitational field at the one- loop order. The structure is a slight generalisation of @, 

[^n-*] (x;x') = L^aplaD *,. + ^^ja 4 9(At 2 - Ar 2 )Q(At) [ln^Ar 2 - At 2 )\ - l] , 

where a and j3 are field dependent constants and D pa and D a ^ a are the derivative operators defined in ^ (if 
a field is massive, then the scalar function on the second line of ([2]) will be different, but on the light cone it will 
reduce to the one given in ^ and in the equation above). The vacuum polarization tensor at higher loops will 
still have the same tensor structure, but scalar functions will generally be different. 

Conversely, the structure of Eq. (TTJ) is such that it is gauge invariant. Since this equation is linear in metric 
perturbations, it also means that it is covariant to the order considered. This can be shown as follows. Firstly, note 
that both the graviton vacuum polarization ([2]) and the Lichnerowicz operator ((3]) consist of operators Z? M „ and D^ vpa 
defined in Eqs. Now, since under an infinitezimal coodinate shift, x^ — > x p +£ M (ai), where £ M (x) is an infinitczimal 
function, to linear order the metric perturbation h pv transforms as h pv — > h pv — 2di^ u y This then implies that, to 
linear order in £^ we have, 

V -> ^ V - ^D^d^ = D^h pv , 

and 

D^hpa -> D^hpa - 2D tlup ' y d p ^ = D^flpa . 

This immediately implies that the classical part of equation (TT]) is gauge invariant. In the quantum part of ([1]) one 
can perform two partial integrations to move D^ 1 " and D pv( " 7 to act on h pa under the x' integral, implying that 
- up to boundary terms - the whole equation is gauge invariant, or covariant up to the order one considers. This 
analysis applies in Minkowski background. In general (curved) backgrounds, coordinate transformations have an 
analogous form, with the difference that ordinary derivatives are replaced by covariant derivatives both in coordinate 
transformations and in the derivative operators in the vacuum polarization. The only subtlety is that one should take 
proper care of operator ordering, since covariant derivatives do not in general commute. Hence, we make the following 
conjecture: 

On general curved spacetime backgrounds the linearized quantum- corrected equation of motion for gravitational 
perturbations is gauge invariant. 

This conjecture should hold for arbitrary quantum fluctuations of matter and interaction fields. Of course, this does 
not say anything about gauge invariance of the vacuum polarization tensor itself. Invariance under some symmetry 
transformation of the vacuum polarization tensor can be checked by considering, for example, the vacuum polarization 
induced by the one-loop (Abelian or non-Abelian) gauge fields, see for example 0, [H| and [27j . Furthermore, it 
would be of interest to investigate whether gauge invariance persists to higher order in metric perturbations, i.e. 
whether the quantum corrected equations for gravitational perturbations are covariant to second or higher order in 
perturbations. 

Finally, it would be of interest to clarify what is the connection (if any) between the results of the asymptotic 
safety program (according to which quantum gravity in the ultraviolet becomes weak, and hence antiscreened) and 
the results presented in this work. 



Appendix: Solving the integral equation for the gravitational field 

Here we show how to solve the integral equation ([6]) for the (quantum) gravitational field h a ^(x), 



30 U 



l\ 2 ' 



G 



where 



Upon performing a spatial Fourier transform, Eq. (|54j) can be reduced to, 



F pa {x) = d 4 J d 4 x'Q(At 2 - Ar 2 )Q(At) [ln(p 2 (At 2 - Ar 2 )) - l] h prT (x'). 



47T 



oo JO 



F pa (k,t) = — (d 2 +k 2 ) 2 / dt' d(Ar)Arsin(fcAr) In U 2 (Ai 2 - Ar 2 ) - 1 h p(T (k, t') , 



F pa {k,t)= d 3 xe- lk - x F pa (x,t), h prT (k,t')= / d 3 xe~ lk x h pr7 {x,t') . 



where k = ||fc|| and 



The spatial integral can be reduced to S-i(z, Q defined in Ref. [261 ]. Eqs. (45-50), 

F pa {k,t) = ^{d 2 + k 2 ) 2 J' dt'E^XKAKt'), 

where z = kAt, £ = k/ fi, x = Ar / At, 

Z n (z,() = z 2 ( dxx sin(za;)[21n(z / 'Q + ln(l - x 2 ) +n]= [sin(z) - z cos(z)] [2 ln(z/Q + n] + z 2 £(z) 
Jo 

and 



£(z) = z 2 / dxa; sin(za;) ln(l — x 2 ) 
Jo 

— 2sin(z) — [cos(z) + z sin(z)] [si(2z) + ir/2] + [sin(z) — zcos(z)][ci(2z) — 7e — ln(z/2)] 



with 



si(z) 



sin(t)dt f z sin(t)dt ir 



t 



t 



ci(z) 



cos(t)dt f z [cos(t) - l]dt 



t 



t 



Now, since for small z = kAt, S n oc z 3 ln(z), one can act with d 2 + k 2 once without touching the upper limit 
time integral in (|56|) to obtain, 

F pff (fc,0 = T (5 t 2 + fc 2 ) / ^'(a 2 + i)s_ 1 (z,c)V(^0, 



where 



(<9 2 + l)S„(z,C) = sin(z)[21n(z/C) + (n + 1) + ci(2z) - lE ~ ln(2z)] - cos(z)[si(2z) + tt/2] 



Next, upon observing that, 



(d 2 + l) 2 E n (z, C) = 4 C ° S(Z) 1 + 4d z [ln(2z/C) + (n + l)/2] , 



one concludes that (|59|) can be recast as 
F pa (k,t) = 16tt 



* tf^^+fl, f dt'ln(2,At) 



h pa (k, t ) 
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We shall now argue that in the second integral, one can move the derivative dt to act as df on h pa (k, t'). Notice firstly 

that one can write ln(2/xAf) = — 9f{At[ln(2//At) — 1]}, and then one can partially integrate df to act on h pr7 (k,t'). 
Formally, this will generate an infinite contribution from the initial state at t' — > — oo. This infinite contribution 



can be regulated away, if one assumes that h pa (k,t') gets adiabatically switched on as h p<7 (k,t') oc e £ * as t' — > — oo. 
After this, one can act with d t , and the second integral becomes f_ dt' \n(2[iAt)dfh pcr (k, t'). This expression is now 
prepared for a Fourier transform over time. By noting that 



•2tt 



dk° 



tk"At 



2tt 



F pa {k p ) can be written in momentum space as, 



FpaW) = -8tt 



In 



|(fc°) 2 -k 2 



2 7£ ; - y (sign(fc° + fc) + sign(fc° - k) 



h p a{k^) ■ 



where we made use of the following integrals, 



\k°f 



At 



sign(k° + k) + sign(/c° - fc) - 2sign(fc°) 



/ dAte lk ° At In(2/iAt) = - — 



l n (^)-^sign(fc°) +7£ 



Inserting (|61j) into Eq. f)53|) and transforming into momentum space one gets, 



h afi{c) {k p ) = h a0 {k p ) + 
x In 

where (cf. Eq. (j4j) 



960tt 2 
|(fc°) 2 -fc 2 | 
4^ 2 



D af3 p ° - Tl a pD pa 



7 + 



1\2 



- 2 7£ ; - y (sign(fc° + fc) + sign(fc° - k) 



h p a(k p ) , 



D p ° = -Ffc ff + TT(k»kJ , D a /" = -k (a 6 \ p k^ + hj p S^(k p k p ) + \^k a k^ 



(61) 



(62) 



(63) 



Next it is useful to act with the two operators (|63|) on h pa - It is interesting to note that these generate only gauge 
invariant combinations of components of h pa . For example, 



D pa h pa (k p ) = 2k 2 <i>(k p ) + 2[3(fc ) 2 - 2k 2 ]^{k p ) , 



(64) 



where 



* = ~(h + k 2 h) , $ = -i(ftoo + 2«fc°(T - {k°) 2 h) . 



and we have made use of a Helmholtz decomposition of h pu , 

5,: 



hoi = n- + ikiO . 



3 



hi kj 



-jtk 2 yh + i(hhj + kjhj) 



(65) 



(66) 



where n\ ,hf and hjj are transverse, k l nj = = k l hj , k l hjj = = hjj k 3 , and /i^- is in addition traceless, 

hfi T = and k = \\k\\. Apart from the Bardeen potentials, there is one gauge invariant vector and one gauge invariant 
tensor, 



n 4 



nj + ik°hj . 



■TT 



(67) 



such that there are in total 6 gauge invariant fields, which means that four out of the ten components of h pv are gauge 
dependent, i.e. they change if one perform a coordinate shift x 1 ' — >• .r** +£, p (x), where £ A '(.?') is a small coordinate 
shift, of the order of metric components h pv . Note that hj? is on its own gauge invariant, and represents the two 
dynamical, tensorial degrees of freedom of gravity (gravitational waves) . 
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Similarly, the action of the second operator in (j63]) yields gauge invariant quantities only, 



D pa h 
D p<T fr 



D pa h 

J^ij 'l>pcr 



-fc 2 $-3(fc u ) 2 $ 



= ~k 2 nf 
2 

= kikj$ - 



+ 2k k i y = D i0 f " T h p 



kj 



5 l3 ( - (fc ) 2 + fc 2 )] * - kPk^hl + i [ - (fc ) 2 + fc 2 ] ft 



TT 



(68) 



Since the Lichnerowicz operator ([3]) of the classical equation of motion for fif) is a linear combination of D p(7 and 
D a p a , it also yields gauge invariant quantities in the equation of motion. This observation was already made in 
our earlier work [l|. There is however a much stronger statement one can make: the structure of the equation that 
includes the graviton vacuum polarization (jTH2j) is also gauge invariant (at least to the order the equation is written) . 
This property can be used to our advantage and study quantum effects in a fully gauge invariant fashion. 

Since all components of the vacuum polarization are gauge invariant, one can use the Helmholtz decomposition (|66[) 
and simply combine the component equations (I62[) into gauge invariant combinations according to (|65l) and (|67l) . 
resulting in the following gauge invariant equations: 



\|/( c ) 



~T(c) 

n. 



(c) _ 



-T 



-ic{(-(k°r 



fc 2 )$ 



fc 2 )*- 

+ £(-(fc°) 2 + fc 2 )nf 



■30(£- 
-30U- 



- 2fc 2 $- 

- 2fc 2 $ 



2(-3(fc ) 2 - 
-2(-3(fc ) 2 



2fc 2 )* |, 
-2fc 2 )* |, 



h TT 



(c) 



hjr + lC(-(k ) 2 + k 2 )h 



■TT 
Hj ■ 



where we defined 



K = 



1920tt 2 



In 



V V 



2 7£ ; - — (sign(fc° + fc) + sign(fc° - fc) 



The last two equations ([711172")) are easily solved, 



-Tfc) 



1 + ( - (fc ) 2 + k 2 )lC(k») 
l + fc(fc")(- (fc°) 2 + fc 2 ) 



(69) 

(70) 
(71) 
(72) 

(73) 

(74) 
(75) 



and tell us that the effect of one loop vacuum polarization for the vector and tensor perturbations can be absorbed 
into a finite wave function renormalization, i.e. vacuum polarization changes the amplitude of the gravitational wave. 
The effect is however tiny, since on shell, where — (fc ) 2 + fc 2 = 0, the effect vanishes. That means that, when there are 
no classical gravitational waves, quantum effects of scalar vacuum fluctuations do not generate any. The same holds 
for vector perturbations: if they are zero classically, they will not be generated by fluctuating scalar fields. Because 
the structure of vacuum polarization is dictated by gauge invariant tensor structures, we expect the same conclusions 
concerning gauge invariant vector and tensor perturbations to hold for arbitrary fields: the quantum effects of the 
fluctuations will just renormalize the wave function. 

The structure of the scalar equations (|SM70"|) is more complex, because they in general couple. Note that the 
equation for the sum of the Bardeen potentials is particularly simple, 



$W + *( c ) = $ + <j + jc( - (fc ) 2 + fc 2 ) ($ + *), 
whose solution is of the same form as that of the vector (|74|) and tensor perturbations l|75[). 

<j>( c ) + \i/( c ) 

$ + * 



1 + /C(- (fc°) 2 + fc 2 ) 



(76) 



(77) 



i.e. the effect of loop fluctuations of quantum fields on the sum $ + ^ is just a finite renormalization of its amplitude, 
which vanishes on-shell. 
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It is, of course, possible to decouple the two Bardeen potentials 
equations is, 



A convenient way of writing the decoupled 



* = 



<j>( c ) 



1 - 5(1 - 6£) 2 (iWb M )/C 

<f( c > 1 2 

1-5(1-60 2 (^M/C ~ 3 fc 



(fc u ) 2 £ 



5(1 ~6£) 2 



1 - 5(1 - Q^f{k^)K 1 + (fc^fcJ/C 
5(1 -6£) 2 1 



1 - 5(1 - 6£) 2 (fc>%)/C l + (fc^fe M )/C 



($(«0 + *( c )) , (78) 
(79) 



This means that each of the two quantum-corrected Bardeed potentials depend on both classical Bardeen potentials. 
Since in general relativity, the classical potentials sourced by a point mass (or a scalar field) are equal $^ c ^ = \l/( c ), 
equations ([75H75)) further simplify to, 



1 1 - 10(1 - 6£) 2 (fc°) 2 /C 2 2 + (fc°) 2 /C 



3 1-5(l-6£) 2 (fc^fc M )/C 3 1 + 

1 2-5(1 -6£) 2 (fc°) 2 /C 1 1 - (fc°) 2 /C 
3 1 - 5(1 - 6£) 2 (fe>%)/C + 3 1 + (/^fc M )/C 



(80) 
(81) 



For a point stationary mass we have $( c ) = = 4nGNM/k 2 . These equations are now used to obtain the static 
limit of a (one-loop) quantum corrected gravitational response to a point mass in the main text. 
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